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Abstract 

An analytic collective model in which the relative presence of the quadrupole and oc- 
tupole deformations is determined by a parameter (0o)> while axial symmetry is obeyed, 
is developed. The model [to be called the analytic quadrupole octupole axially symmet- 
ric model (AQOA)] involves an infinite well potential, provides predictions for energy and 
B(EL) ratios which depend only on (p , draws the border between the regions of octupole 
deformation and octupole vibrations in an essentially parameter-independent way, and de- 
scribes well 226 Th and 226 Ra, for which experimental energy data are shown to suggest 
that they lie close to this border. The similarity of the AQOA results with <p = 45° for 
ground state band spectra and B(E2) transition rates to the predictions of the X(5) model 
is pointed out. Analytic solutions are also obtained for Davidson potentials of the form 
(3 2 + Pq/P 2 , leading to the AQOA spectrum through a variational procedure. 
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1. Introduction 

Rotational nuclear spectra have long been attributed to quadrupole deformations pQ, 
while octupole deformations [corresponding to reflection asymmetric (pearlike) shapes] are 
supposed to occur in certain regions, most notably in the light actinides [21 OH EJ [H] . The 
hallmark of octupole deformation is a negative parity band with levels L n = 1~, 3~, 5~, . . . , 
lying close to the ground state band and forming with it a single band with L 71 " = + , 1~, 
2 + , 3~, 4 + , . . . , while a negative parity band lying systematically higher than the ground 
state band is a footprint of octupole vibrations. The transition from the regime of octupole 
vibrations into the region of octupole deformation has been considered by several authors 
(HI IE]. A complete algebraic classification of the states occuring in the simultaneous 
presence of the quadrupole and octupole degrees of freedom has been provided in terms of 
the spdf- interacting boson model [HI GDI, involving free parameters. (It should be noted, 
however, that an alternative interpretation of the low-lying negative parity states in the 
light actinides has been provided in terms of clustering [HI [T2J US]-) 

On the other hand, the transitions from vibrational [U(5)] shapes to axially symmetric 
deformed [SU(3)] and 7-unstable deformed [SO(6)] shapes have been recently described 
in terms of the X(5) and E(5) ^2] models respectively, which utilize an infinite well 
potential in the (3 degree of freedom, leading to parameter-free (up to overall scale factors) 
predictions for spectra and transition probabilities. 

It is the aim of the present work to provide an analytic description of the light actinides 
lying near the border between the regions of octupole vibrations and octupole deformation, 
through the use of a model containing the minimum number of free parameters. In this 
direction, the following steps are taken: 

1) Quadrupole and octupole deformations are taken into account on equal footing, their 
relative presence decided by the only free parameter in the model, (f)Q. 

2) Axial symmetry is assumed, in order to keep the problem tractable. 

3) Symmetrization of the wave functions is carried out as in Ref. [ISj, involving the 
irreducible representation (irrep) A of the group D 2 for the levels of even parity and the 
irrep B± of the same group for the levels of odd parity. 

4) Separation of variables is achieved in a way analogous to the one used in the framework 
of the X(5) model [Hj. 

5) An infinite well potential is assumed appropriate for the description of the border 
region, as in the E(5) ^Hj and X(5) ^1] models. 

The predictions of the model, to be called AQOA, are compared to spectra and B(EL) 
ratios for 226 Th and 226 Ra, for which evidence from systematics of experimental data is 
presented, suggesting that they lie close to the border between octupole deformation and 
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octupole vibrations. This border is found to be drawn by the AQOA model in an essentially 
parameter-independent way. 

In addition, solutions for Davidson potentials |T7j of the form (3 2 + Pq/P 2 are obtained, 
and a variational method [TSJ E]] leading from the Davidson results to the AQOA predic- 
tions is worked out. 

A different approach to the problem of phase transition in the octupole mode has been 
recently given in Ref . [20] , where the starting point is the introduction of a new parametriza- 
tion of the quadrupole and octupole degrees of freedom, using as intrinsic frame of reference 
the principal axes of the overall tensor of inertia, as resulting from the combined quadrupole 
and octupole deformation. Comparisons between the results of the two methods are de- 
ferred to the appropriate sections. Three main differences between the two models are: 

1) The AQOA model is analytic, while the model of Ref. is not. 

2) In the AQOA model the quadrupole and octupole degrees of freedom are taken 
into account on equal footing, while in the special form of the model of Ref. [20] used 
for comparison to experiment, the octupole degree of freedom remains active, while the 
quadrupole degree of freedom is "frozen" to a constant value. 

3) In the AQOA model the symmetry axes of the quadrupole and octupole deformations 
are taken to coincide, in order to guarantee axial symmetry, while in the more general 
framework of Ref. [20] nonaxial contributions, small but not frozen to zero, are taken into 
account. 

In Section 2 the AQOA model is formulated, while numerical results are given in Section 
3 and compared to experiment in Section 4. In Section 5 the variational procedure is 
described, while Section 6 contains discussion of the present results and plans for further 
work. 

2. The Analytic Quadrupole Octupole Axially Symmetric (AQOA) Model 
2.1 Formulation 

We consider a nucleus in which quadrupole deformation (flz) an d octupole deformation 
(^3) coexist. We take only axially symmetric deformations into account, which implies that 
the 7 degrees of freedom are ignored, as in the Davydov-Chaban approach j2I] . The body- 
fixed axes x', y', z' are taken along the principal axes of inertia of the (axially symmetric) 
nucleus, while their orientation relative to the laboratory-fixed axes x, y, z is described by 
the Euler angles 9 = {9\, 92, #3}. The Hamiltonian reads ^01 122j 



E 



h 2 1 d 3 d 



h 2 L 2 





H 



2B x (3ldf3 x HX d(3 x 



+ 



6(B 2 (3 2 + 2B 3 (3 2 ) 



A=2,3 



where B 2} B 3 are the mass parameters. 
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We seek solutions of the Schrodinger equation of the form jTH] 

QHfaM = (/3 2 /3 3 )- 3/2 vI/±(/5 2 ,/33)|LM0,±), (2) 

where the function |LM0,±) describes the rotation of an axially symmetric nucleus with 
angular momentum projection M onto the laboratory-fixed z-axis and projection K = 
onto the body-fixed z'-axis. The moment of inertia with respect to the symmetry axis z' 
is zero, implying that levels with K 7^ lie infinitely high in energy [16J. Therefore in 
this model we are restricted to states with K = only. The function \LM0, +) transforms 
according to the irreducible representation (irrep) A of the group D 2 , while the function 
\LM0, — ) transforms according to the irrep Bi of the same group fnj|22]- The general 
form of these functions is £Q 



\LMK, ±) = J_^±L_(2^ iAr (e) ± (-1) L V L KM (6)). (3) 



In the special case of K = it is clear that \LM0, +) 7^ for L = 0, 2, 4, . . . , while 
\LM0, — ) 7^ for L = 1, 3, 5, . . . The functions ^^(fo, P3) and ^2(^2,^3) are respectively 
symmetric and antisymmetric with respect to reflection in the plane x'y', and therefore 
describe states with positive and negative parity respectively [2~2~] . 

Using the solutions of Eq. (J2J) for the Hamiltonian of Eq. ((TJ) the Schrodinger equation 
takes the simplified form 

h 2 d 2 h 2 d 2 h 2 L(L + i) 



2B 2 d(3 2 2B 3 d(3 2 6(B 2 PI + 2B 3 P!) 



3h 2 



This equation is further simplified by introducing [TBJ 122] 



=0. (4) 




$2=fi2\\— , /3a = /3sy B= , (5) 

as well as reduced energies e = (2B/K 2 )E and reduced potentials u = (2B/h 2 )V [11113, 
reaching the form 



d 2 d 2 L(L + 1) - ~ 3/1 1 



d Ji-Wl + W^) +u{ ^ ) + ~AM + W~ eh 



= 0. (6) 



Further simplification occurs through the introduction of polar coordinates (with < (3 < 
00 and -tt/2 < < tt/2) [E1I22] 



/5 2 = /?cos0, /3 3 = /?sin0, (3 = ^[3$ + (3l (7) 
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leading to 

dp 2 



1 d L(L + 1) 1 d 2 

pd/3 + 3/3 2 (l + sin 2 0) ~ /f 2 d^ 2 



+ «(/3,0) + -~ 



f3 2 sin 2 20 



It is clear that = corresponds to quadrupole deformation alone, while = ±tt/2 
corresponds to octupole deformation alone. It is worth noticing that the transformation of 
Eq. (J7J) allows (3 3 to assume both positive and negative values, while (3 2 takes only positive 
values. 

Separation of variables in Eq. (jHJ) can be achieved by assuming the potential to be of 
the form u(/3, 0) = u{(3) + m(0 ± ), where m(0 ± ) is supposed to be of the form of two very 
steep harmonic oscillators centered at the values ±0o, i.e. 



u(4> 



T0oj 



r±\2 



=F 0o, 



(9) 



with c being a large constant. In other words, the nucleus is supposed to be rigid with 
respect to the variable 0, implying that remains close to ±0o and, therefore, the relative 
amount of quadrupole and octupole deformation remains constant, as in Strutinsky-type 
potential energy calculations 0. This assumption will be (partly) justified a posteriori by 
the fact that the spectrum remains almost unchanged for values of 0o between 30° and 60°. 
In this way Eq. (jHJ) is separated into 



d 2 



d(3 2 



1 d 1 

Pdp + J 2 



' L(L + 1) 
3(1 + sin 2 0, 



+ 



sin 20 o 



+ U 0)-6g(L) 



= 0, (10) 



and 



d 2 



+ u( 



- e„ 



x{4> 



0. 



((3 2 ) 9(0±) 2 

where *t(A0) = ipt(P)(x(<t> + ) ± x(0"))/v / 2, while {(3 2 ) is the average of (3 2 over ^0), 
and €l = e p{L) + e</,. It is worth noticing that Eq. (fTU|) has the same form for both +0o 
and —0o, since only even functions of 0o appear in it. 
2.2 The /3-part of the spectrum 

In the case in which u([3) is an infinite well potential {u{f3) — if (3 < $w; u{(3) = oo 
if (3 > Pw)i using the definitions = k\, z — $kp, Eq. (fTUj) is brought into the form of a 
Bessel equation 

d 2 ^ , 1#± 



c^ 2 



+ 



+ 



with 



L(L + l) 



\ 3(1 + 



sm 



+ 



sin 20o 



(12) 



(13) 
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Then the boundary condition ip^(Pw) — determines the spectrum 



e (3,s,v e f3,s,<f> ,L 



) 2 
,v) ) 



k 







(14) 



w 



(15) 



and the eigenfunctions 

where x SyV is the sth zero of the Bessel function J v {z), while c SjU are normalization constants, 
determined from the condition Jq \ipf^((3)\ 2 (3d(3 = 1 to be c SjU = \/2/ J u +i(k s ^). The 
notation has been kept similar to Ref. |14j . 

Eq. (fTUj) is also exactly soluble (23 121] in the case of the Davidson potentials [Uj 



(16) 



In this case the second term of Eq. (fTH|) is combined with the third term of Eq. (flOjl . 
leading to eigenfunctions which are Laguerre polynomials 



2n\ 



T(n + a+ 1) 



(17) 



where 



L(L + 1) 



+ 



^3(l + sin 2 o ) sin 2 20o 
while the energy eigenvalues are given by 



E nL = 2n + a + l = 2n + l + 



L(L + 1) 



\ 3(1 + 



sin O ) sin 20, 



(19) 



It is worth remarking that the excitation energies, Eq^—Eqq, within the ground state band 
(which is characterized by n = 0), divided by an appropriate normalization constant read 



E' 



0,L,exc 



1 + b H L(L + 1) - 1, 



(20) 



with b H l = 3(1 + sin 2 O ) ( si J 2(j>0 + Eq. fl2D|) is the Holmberg-Lipas formula |2S1- 
In what follows, the infinite well potential will be used everywhere. Davidson potentials 

will be briefly employed in Section 5. 
2.3 The 0-part of the spectrum 
Eq. ()11|) for the potential of Eq. Q takes the form 

d 2 



(21) 
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where =t = <p =p 0o- This is a simple harmonic oscillator equation with energy eigenvalues 
and eigenfunctions 



2c / 1\ 

n<t> + ~ ), n = 0,1,2,... (22) 



\ (P 2 ) 



2 



Y,#) = iV^^e"^ 2 / 2 , 6 = ('^V /4 , (23) 



with normalization constant AL, = . / 7 =-4i — r . 
The total energy in the present model is then 

E(s, L, O , n<p ) = E + Ae^ L + Bn^. (24) 

2.4 B(EL) transition rates 

In the axial case used here the electric quadrupole and octupole operators are 

Tf 2) = t 2 [3 2 V^(9), Tf 3 ) = t 3 p 3 vfl(9), (25) 
while the electric dipole operator reads [TB] 

T^^ttMP^iO). (26) 
The total wave function in the case of the infinite well potential is 



(27) 

where C is a constant, while in the case of the Davidson potentials the same expression 
holds with J u (k s ^ u f3) replaced by F£(/3). 
B(EL) transition rates are given by 

8[H;te , V)) . i|3M , (28) 

where the reduced matrix element is obtained through the Wigner-Eckart theorem 

(L fm \T^\L^) = (Wfl^f) {L Q/ || T ^|| W (29) 

/2L f + 1 



In Eq. f!28j) the integration over the angles # involves a standard integral over three 
Wigner functions [26J, which leads to (LjLL^|000), while the rest of the integrations are 
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performed over / / /Jfd/^/Jfci/^, where the /3f factors come from the volume element 
and cancel with the first factor of Eq. (J2*7|) . Using Eqs. (0) and (J2J), as well as the relevant 
Jacobian, one finds (up to constant factors) that the integration is over / (3df3d(f). 

In the integrals over 0, only the case of = 0, corresponding to H = 1, is considered. 
The results are factors depending on the parameters b and </>o, as well as on the multipolarity 
of the transition. Therefore in Section 3, ratios of B(EL) transition rates will be presented, 
in which these factors cancel out. 

The integrals over /3 are 

if 2) = if 3) = J ^J Ui (k Si ,J)J Vf (k Sf>Vf ^, (30) 

lf l) = J p 3 J Ui (k Sl ,J)J Uf (k Sf , Uf P)dp, (31) 

in the case of the infinite well potential, while for the Davidson potentials the Bessel func- 
tions are replaced by Laguerre polynomials, as above. The final result then reads 

B(EL; U^L S ) = c(c s ^ i c S/ ,, / ) 2 (L J LL / |000) 2 (/^ L) ) 2 , (32) 

where L = 1,2,3 and all constant factors have been absorbed in c. 
3. Numerical results 

Spectra for the ground state band and the negative parity band associated with it 
(s — 1), as well as for the first excited band (s = 2) and the second excited band (s = 3), 
normalized to the 2f state of the ground state band, are shown for several values of 0o 
in Table 1. A few R(L) = E(L)/E(2) ratios are also depicted as functions of 0o i n Fig. 
1(a). It is clear that the results are quite stable in the region 30° < 0o < 60°, while at the 
limiting cases near O — 0° and 90° the rigid rotor results are obtained, corresponding to 
a pure rotational spectrum for the ground state band and the associated negative parity 
band, while the excited bands are pushed to infinity. 

B(E2) transition rates are listed in Table 2 for several values of <fio, while a few B(E2) 
ratios are shown in Fig. 1(b) as functions of <fio, their behavior being quite smooth in the 
region 30° < 0o < 60°. The same remark applies to B(E1) and B(E3) transitions, listed 
in Tables 3 and 4, and shown in Fig. 1(c). 

It is worth remarking that the minima of energy ratios related to the ground state band, 
as well as the maxima of B(EL) ratios regarding the ground state band and the associated 
negative parity band, reported in the caption of Fig. 1, are all located between 0o — 40° 
and 43°, while the minima of energy ratios regarding the excited (s = 2, 3) bands are 
located near O — 35°. 
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The tables and figures mentioned so far indicate that the region of interest in the present 
model, in which smooth and essentially parameter independent behavior of spectra and 
B(EL) rates is observed, is the region 30° < 4>q < 60°, to which further considerations will 
be limited. 

In addition to the results of the AQOA model, the X(5) spectrum is included in Table 
1 for comparison. It is clear that the ground state band of X(5) lies a little lower than 
the ground state band of the AQOA model with 0o — 45°, while for the s = 2 and s = 3 
bands the AQOA model predictions for 0o = 45° are larger than the X(5) values by almost 
a factor of two. Furthermore, in Table 2 the B(E2) transitions within the ground state 
band of X(5) are shown for comparison. It is clear that the X(5) values are slightly higher 
than the corresponding predictions of the AQOA model for <p = 45°. 

The similarities between the ground state bands of the AQOA and X(5) models can be 
understood as due to the fact that both models originate from the Bohr Hamiltonian and 
use an infinite well potential, while in addition for the properties of the ground state band 
the quadrupole degree of freedom, included in both models, is expected to be important. 
In contrast, the excited bands appear to be more sensitive to the inclusion of the octupole 
degree of freedom. The position of the state becomes therefore an important factor in 
the process of comparison to experiment. One can also think of the AQOA model as an 
extension of the X(5) framework, in which the negative parity states, as well as the B(EL) 
transitions involving them, are included. 

4. Comparison to experiment 

Experimental data for the ground state and related negative parity bands of 220_234 Th 
are shown in Fig. 2(a). It is clear that 226 Th lies on the border between two different 
regions. Below 226 Th the odd-even staggering is very small, while from 228 Th up the 
odd-even staggering is becoming much larger, increasing with the neutron number N. A 
quantitative measure of the odd-even staggering and related figures can be found in Ref. 
[37J. It is clear that below 226 Th the situation corresponds to octupole deformation, in 
which the ground state band and the negative parity band merge into a single band, while 
above 226 Th the picture is corresponding to octupole vibrations, i.e. the negative parity 
band is a rotational band built on an octupole bandhead, thus lying systematically higher 
than the ground state band. Theoretical predictions for = 45° lie a little below 226 Th, 
while the <fi = 60° results follow the 226 Th data very closely. It is worth remarking that 
the procedure of Ref. [20 , which is quite different from the present one, also leads to the 
identification of 226 Th as the nucleus lying closest to the transition point from octupole 
deformation to octupole vibrations. 
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A similar picture is observed in 218_228 Ra, shown in Fig. 2(b). In this case octupole 
deformation appears below 226 Ra, while 228 Ra is already in the regime of octupole vibrations. 
Theoretical predictions for 0o — 45° again lie a little below 226 Ra, while the 226 Ra data are 
followed quite closely by the predictions of 0o = 56°. 

The behavior observed in Fig. 2(a) can be better understood by considering Figs. 3(a) 
and 3(b), where the experimental energy levels of the ground state band and the associated 
octupole band are shown, for the same thorium isotopes. While the even parity levels, shown 
in Fig. 3(a), smoothly decrease with increasing neutron number N, as a result of increasing 
quadrupole collectivity, the odd parity levels, shown in Fig. 3(b), exhibit a minimum, which 
is located at iV = 136 up to L = 9, while it moves to N = 138 for higher L. This change 
of behavior is then attributed to the octupole degree of freedom, showing that gQ 6 Thi 36 
lies near the border between octupole deformation and octupole vibrations. The change of 
behavior is not abrupt, since the effect due to octupole deformation is "moderated" by the 
quadrupole deformation setting in in parallel. 

In a similar manner the behavior observed in Fig. 2(b) can be clarified by consider- 
ing Figs. 3(c) and 3(d), where the experimental data for the same radium isotopes are 
presented. Again, the even parity levels decrease with increasing N, while the odd parity 
levels exhibit a minimum, located at iV = 136 up to L = 5, while it moves to iV = 138 
for higher L, showing that gg 6 Rai 38 lies close to the border between the regions of octupole 
deformation and octupole vibrations. 

The transition from octupole deformation to octupole vibrations can also be seen by 
considering the simplest quantity measuring the relative displacement of the negative parity 
levels with respect to the even parity ones, 

AEiL) = E(L) - E(L-l) + E(L + l) (33) 

Results for the Th and Ra isotopes are shown in Figs. 4(a) and 4(b), respectively. In Fig. 
4(a) it is clear that in 222 - 22 6Th the stag gering is decreasing rapidly with increasing angular 
momentum, reaching a vanishing value and staying close to it, which is the hallmark of 
octupole deformation [21EH1, while in 228 ~ 234 Th the decrease is much slower and vanishing 
values, if any, correspond to very high angular momenta, a behavior expected for octupole 
vibrations. Again 226 Th appears closest to the border between the two regions. In Fig. 
4(b), 220-226 Ra exhibit the rapid decrease of staggering and the sticking to values close to 
zero beyond the first vanishing value, while 228_230 Ra follow the slow decrease pattern. As 
a result, 226 Ra appears to be closest to the border line between the two regions. 

As far as the bandhead is concerned, the experimental values (normalized to the 2f 
state) are 12.186 for 226 Ra and 11.152 for 226 Th, in good agreement with the 11.226 and 
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12.410 values predicted by the AQOA model for the 0o values of 56° and 60° used in Fig. 
2 . (The model of Ref. [20] provides a value of 8.528 for 226 Th.) It should be noticed 
that the normalized bandhead is lying close to this height for all Ra and Th isotopes 
for which data exist, namely 222 Ra (8.225), 224 Ra (10.861), 228 Ra (11.300), 228 Th (14.402), 
230 Th (11.934), 232 Th (14.794), 234 Th (16.347), with data taken from the references used in 
Fig. 2 . 

Considering the AQOA model as an extension of the X(5) framework involving negative 
parity states, as remarked at the end of Section 3, implies that the search for X(5)-like nuclei 
in the light actinides, where the presence of low- lying negative parity bands is important, 
should be focused on nuclei with i?(4) ratio close to 3.0 and 0^ bandhead higher than the 
X(5) value of 5.65 . 

Detailed comparisons to BiEV) transition rates are not feasible, because of lack of 
experimental data. We therefore use ratios of B(EL) transitions, also used in earlier work 
[201 ESI- Thus in Table 5 and Fi S- 5 ( a ) the experimental B(E1;L — > L — 1)/B(E2;L 
L — 2) ratios used in Ref. [20] are shown, together with theoretical predictions from the same 
source, and predictions for <p = 45° and 60°, the values also used in Fig. 2(a). The present 
theoretical predictions for the two different values of 0o practically coincide (indicating that 
the predictions are essentially parameter free) and are in most cases within the error bars of 
the experimental points, while the predictions of Ref. [20] grow a little faster as a function 
of angular momentum. 

Furthermore, in Table 6 and Fig. 5(b) the experimental B(E1\ L —>■ L + 1)/B(E1; L — ► 
L — 1) ratios J30] used in Ref. [2H| are shown, together with three sets of theoretical 
predictions in the framework of the extended coherent states model (ECSM) jHj from the 
same source, corresponding to the lowest order choice for the El transition operator (R- 
h), as well as to two different choices of the El transition operator, including anharmonic 
terms assumed suitable for the transition region (R-I, R-II) [SHI- In addition, predictions 
for 0o — 45° and 56°, the same values used in Fig. 2(b), are shown. It is clear that 
the predictions for the two different values of 0o practically coincide (indicating that the 
predictions are essentially parameter free) and in all cases are within the error bars of the 
experimental points, being in very close agreement to the R-I predictions of Ref. [3*9"j . 

On the results presented in this section, the following additional comments apply. 

1) Figs. 2 and 5 indicate that 226 Th ( 226 Ra) can be well described using the AQOA 
model with (f> = 60° (0o = 56°), which provides results quite similar to the 0o = 45° case. 
In all these cases, Eq. (J7J) [together with Eq. (jSJ)] indicates that the quadrupole and octupole 
deformations are present in comparable amounts. This is in agreement with Strutinsky- 
type potential-energy calculations 01121, resulting in comparable @ 2 an d 03 values for these 
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nuclei. The presence of octupole deformation in 226 Ra has also been realized in a study jlHI 
within the framework of the spdf-IBM |H| ll()j. 

2) Figs. 2-4 suggest that 226 Th and 226 Ra lie close to the border between octupole 
deformation and octupole vibrations. This is in agreement with Woods-Saxon-Bogolyubov 
cranking calculations [7j for the Ra and Th isotopes, suggesting shape changes from nearly 
spherical (N ~ 130) to octupole-deformed (N ~ 134) to well-deformed reflection-symmetric 
(N ~ 140) shapes, in which negative-parity bands can be interpreted in terms of octupole 
vibrations. 

3) One can easily see that no odd-even staggering is predicted by the AQOA model. 
This is in agreement to the well known fact that odd-even staggering is produced when 
the potential in /3 3 is a double well with two symmetric minima |44j . the staggering being 
sensitive to the angular momentum dependence of the height of the potential barrier [4*3] . 
An infinitely high barrier leads to no odd-even staggering |44j . which is indeed the case 
here. The introduction of a finite barrier in the present model will lead to staggering, but 
it will require the addition of at least one new parameter, in contrast to the main goal 
of the present work, which is the description of the border between octupole deformation 
and octupole vibrations with the minimum number of parameters possible. As shown in 
Figs. 2(a) and 2(b), the model does predict the border between the regions of octupole 
deformation and octupole vibrations in an essentially parameter independent way. 

4) It should be noticed that the transition examined here is the one from octupole 
deformation to octupole vibrations as a function of the neutron number in a chain of 
isotopes, which is different from the gradual setting in of octupole deformation as a function 
of angular momentum in a given nucleus, usually studied by considering the odd-even 
staggering EH] , as already discussed in relation to Fig. 4. 

5. The variational procedure 

In Refs. jTHl EI] a variational procedure has been introduced, leading from the results 
of one-parameter Davidson potentials to the parameter-free E(5) and X(5) predictions. 
The same procedure can be applied in the present case, by considering (for given </> ) 
the R(L) = E(L)/E(2) ratios predicted by the Davidson potentials of Eq. (JTfij) for the 
excitation energies of the ground state band and the associated negative parity band, and 
determining for each value of L separately the value of the parameter /3q at which the 
derivative of the ratio R(L) with respect to (3$ has a sharp maximum. The collection of 
R(L) values selected in this way (for the case of <ft — 45°) is shown in Table 7 and Fig. 6, 
together with the limiting cases of (3q = (a vibrator) and (3q — > oo (a rigid rotor). It is 
clear that the collection created through the variational procedure practicaly coincides with 
the predictions of the present model utilizing an infinite well potential, thus indicating that 
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the choice of the infinite well potential indeed correponds to the transition point between 
a vibrator (/3 = 0) and a rigid rotor (j3 — > oo), since it is at the transition point that the 
rate of change of the R(L) ratios is expected to become maximum. 
6. Discussion 

The analytic quadrupole octupole axially symmetric (AQOA) model introduced in this 
work describes well the border between octupole deformation and octupole vibrations in 
the light actinides, which corresponds to 226 Th and 226 Ra in the Th and Ra isotopic chains 
respectively. Some of the main ingredients of the present model, such as the infinite well 
potential and the approximate separation of variables, strongly resemble the ones used in 
the X(5) model, describing the critical point of the shape phase transition from vibrational 
to axially deformed rotational nuclei |14j . determined through the study of potential energy 
surfaces derived from the Hamiltonian of the Interacting Boson Model [IB]. An interesting 
task is the study of the potential energy surfaces resulting in the spdf-IBM jHl HHj, the 
version of IBM including the octupole degree of freedom in addition to the quadrupole one, 
which can possibly lead to the determination of a shape phase transition from octupole 
deformation to octupole vibrations, in a manner similar to the determination of the critical 
point between the spherical and triaxial shapes found recently through the study of the 
potential energy surfaces resulting from an IBM-2 Hamiltonian |3B|- Although some 
early results are given in Ref. [10j . this task is far from complete. The persistence of axial 
symmetry, as well as the importance of parity projection in this context have been empha- 
sized |49| I5U|. The inclusion of staggering in the present model, as well as its application to 
the rare earth region near A = 150, where octupole deformation is known to occur [3J 
are also of interest. 
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Table 1: Spectra of the AQOA model for the ground state band and the associated negative 
parity band (s = 1), as well as for the first excited band (s = 2) and the second excited band 
(s = 3), normalized to the energy of the 2f state, for different values of <fi . The second 
column contains the values obtained slightly above 0° or slightly below 90°. In addition, 
the X(5) spectrum is shown for comparison. See Section 3 for further discussion. 

0o 0°, 90° 15° 30° 45° 60° 75° X(5) 

1/ 

s = 1 

0+ 0.000 0.000 

1- 0.333 0.337 

2+ 1.000 1.000 

3- 2.000 1.969 

4+ 3.333 3.221 

5" 5.000 4.734 

6+ 7.000 6.490 

7- 9.333 8.471 

8+ 12.000 10.666 

9- 15.000 13.065 

10+ 18.333 15.659 

11- 22.000 18.443 

12+ 26.000 21.410 

13- 30.333 24.557 

14+ 35.000 27.881 

15- 40.000 31.379 

16+ 45.333 35.048 

17- 51.000 38.886 

18+ 57.000 42.892 

19- 63.333 47.064 

20+ 70.000 51.402 



s = 2 














0+ 


13.292 


8.983 


9.351 


12.410 


23.896 


5.649 


2+ 


14.893 


10.726 


11.133 


14.160 


25.502 


7.450 


4+ 


18.384 


14.204 


14.630 


17.763 


29.098 


10.689 


6+ 


23.392 


18.820 


19.209 


22.649 


34.410 


14.751 


s = 3 
0+ 


30.940 


21.944 


23.114 


30.316 


55.625 


14.119 



0.000 
0.344 
1.000 
1.921 
3.069 
4.414 
5.935 
7.620 
9.459 
11.445 
13.574 
15.841 
18.245 
20.782 
23.452 
26.251 
29.180 
32.237 
35.421 
38.731 
42.166 



0.000 
0.346 
1.000 
1.912 
3.039 
4.351 
5.829 
7.459 
9.233 
11.144 
13.187 
15.359 
17.658 
20.081 
22.626 
25.293 
28.080 
30.985 
34.009 
37.150 
40.408 



0.000 
0.342 
1.000 
1.938 
3.119 
4.513 
6.098 
7.857 
9.779 
11.857 
14.082 
16.451 
18.959 
21.605 
24.384 
27.297 
30.340 
33.513 
36.814 
40.242 
43.796 



0.000 
0.336 
1.000 
1.981 
3.264 
4.832 
6.667 
8.755 
11.082 
13.635 
16.406 
19.386 
22.567 
25.943 
29.510 
33.264 
37.200 
41.315 
45.607 
50.074 
54.713 



0.000 
1.000 
2.904 
5.430 
8.483 
12.027 
16.041 
20.514 
25.437 
30.804 
36.611 
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Table 2: B(E2; Li — > Lf) values between states of the AQOA model with s — 1. B(E2)s 
with L, and Lf even are normalized to the 2f — > + transition, while B(E2)s with Lj and 
odd are normalized to the 3^f — > l^f transition. The X(5) results are also shown for 
comparison. See Section 3 for further discussion. 







1 o 
1 


1 K° 
10 


oU 


A K° 
40 


OU 


1 




LI 


T 7T 

L f 
















2+ 


0+ 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


4+ 


2+ 


1.429 


1.475 


1.530 


1.539 


1.509 


1.457 


1.599 


6+ 


4+ 


1.574 


1.701 


1.834 


1.855 


1.786 


1.656 


1.982 


8+ 


6+ 


1.648 


1.870 


2.072 


2.104 


2.005 


1.797 


2.276 


10+ 


8+ 


1.693 


2.011 


2.268 


2.309 


2.187 


1.915 


2.509 


12+ 


10+ 


1.723 


2.131 


2.431 


2.480 


2.342 


2.019 


2.697 


14+ 


12+ 


1.746 


2.236 


2.569 


2.626 


2.476 


2.111 


2.854 


16+ 


14+ 


1.762 


2.327 


2.687 


2.751 


2.593 


2.194 


2.987 


18+ 


16+ 


1.776 


2.407 


2.790 


2.860 


2.695 


2.269 


3.101 


20+ 


18+ 


1.787 


2.478 


2.881 


2.955 


2.785 


2.337 


3.200 


3" 


\- 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 




5~ 


3" 


1.179 


1.227 


1.278 


1.286 


1.260 


1.210 




7" 


5~ 


1.257 


1.374 


1.479 


1.496 


1.444 


1.335 




9~ 


7" 


1.301 


1.491 


1.642 


1.665 


1.595 


1.433 




11- 


9~ 


1.329 


1.591 


1.777 


1.806 


1.723 


1.518 




13" 


11" 


1.350 


1.677 


1.890 


1.924 


1.832 


1.594 




15" 


13" 


1.365 


1.752 


1.986 


2.026 


1.927 


1.661 




17" 


15" 


1.376 


1.817 


2.069 


2.114 


2.010 


1.722 




19- 


17- 


1.386 


1.875 


2.142 


2.190 


2.082 


1.777 
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Table 3: B(El;Li — > Lf) values between states with s = 1, normalized to the 1 
transition. 





00 




1° 




15° 




30° 




45° 




60° 




75° 


\- 


0+ 


1 


.000 


1 


.000 


1 


.000 


1 


.000 


1 


.000 


1 


.000 


2+ 


1- 


1 


.200 


1 


.227 


1 


.264 


1 


.269 


1 


.247 


1 


.215 


3~ 


2+ 


1 


.286 


1 


.358 


1 


.455 


1 


.469 


1 


.414 


1 


.328 


4+ 


3~ 


1 


.334 


1 


.467 


1 


.633 


1 


.657 


1 


.564 


1 


.413 


5~ 


4+ 


1 


.364 


1 


.570 


1 


.807 


1 


.842 


1 


.712 


1 


.489 


6+ 


5~ 


1 


.385 


1 


.670 


1 


.977 


2. 


.022 


1 


.857 


1 


.562 


7" 


6+ 


1 


.401 


1 


.768 


2. 


.141 


2 


,196 


2, 


,000 


1 


.634 


8+ 
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1 


.413 


1 


.864 


2 


.299 


2. 


,364 


2, 


,139 


1 


.705 


9" 


8+ 


1 
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1 


.957 


2. 


.449 


2. 


,523 


2, 


,273 


1 
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2 
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2 


.592 


2. 
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.403 


1 


.847 
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.917 
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.334 
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2 
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.411 


3 
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.171 
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3 
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.360 
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.598 


3 


,757 


3 


,351 


2 


.417 


20+ 


19" 


1 


.471 


2 


.778 


3 


.685 


3 


,852 


3 


,436 


2 


.471 
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Table 4: B(E3;Li — > Lf) values between states with s = 
transition. 



1, normalized to the 3 









1 ° 
1 


10 


Qn° 
ou 




/i 

40 




ou 




7p;o 
( 


TTC 
L i 


TTC 
L ,f 






















Q — 
O 


U 


i 


nnn 

.uuu 


1 nnn 
l.UUU 


i nnn 
l.UUU 


i 


nnn 
.UUU 


i 


nnn 
UUU 


i 


nnn 
.UUU 


A + 
4 


1 


i 


QQQ 
.OOO 


1 QC1 

l.oOl 


1 QfiO 

i.ooy 


i 


Q7Q 

.6(6 


i 


,o04 


i 


QAK 

,o40 


cr — 



9 + 
Z 


i 


.010 


l.Ooo 


l.olo 


i 


fi9Q 

.0zo 


i 


,oyo 


i 


KA 7 

.04 ( 





Q — 
O 


1 


.OOZ 


1 71 n 

i. ( iy 


1 ann 

i.ouy 


1 


oor 
.OZO 


1 


77Q 
. 1 IO 


1 


,oyu 


( 


A + 
4 


i 


71 /I 
. ( 14 


1 Q/l 7 

l.o4 ( 


1 070 

i.y ; y 


o 
Z 


nni 
.UU1 


i 


QQ /I 

.yo4 


i 


an/i 
,oU4 


8 + 

o 


cr — 



i 


77/1 

. / /4 


i.yoo 


9 i Qn 
z.loU 


o 
Z 


.loo 


o 

z 


Q7Q 

.U 16 


i 


nnn 

,yuu 


o— 

y 





i 


891 
.oZl 


9 n^s 
z.UOo 


9 9fi7 
Z.ZO / 


o 
Z 


Qnn 
.oUU 


o 

z 


1 no 

. iyy 


i 


oaf; 

.yoo 


i n+ 
1U 


( 


i 


.soy 


9 1 /IO 

z.i4y 


9 Qon 

z.oyu 


2 


.429 


2. 


.314 


o 

z 


nfiQ 
,U0o 


i i - 
1 1 


o 


i 
i 


88Q 
.coy 


9 933 


9 ^n9 


2 


.546 


2 


,420 


9 


1 Qc; 
. lOO 


1 9+ 
1Z 


Q~ 

y 


i 


.y io 


9 qi n 

Z.olU 


9 fine; 
z.ouo 


2 


.653 


2 


,517 


Z 


on9 
,zuz 


lo 


1U 


i 


.yoo 


9 Q81 
Z.ool 


9 P.QQ 

z.oyy 


2 


.752 


2. 


.607 


Z 


OP.A 
,Z04 


1 A+ 
14 


1 1 - 
1 1 


i 


.yoo 


9 AA7 
Z.44 1 


9 78fi 
Z. I oO 


2 


.842 


2. 


.691 


9 
Z 


Q9Q 
OZo 


10 


1 9+ 
1Z 


i 


.y 1 1 


9 c;nQ 
z.ouy 


Z.oOO 


2 


.926 


2 


,768 


9 
Z 


Q7Q 

o i y 


ID 


1 Q~ 
lo 


i 


.yoo 


Z.00 ( 


9 QQQ 

z.yoy 


3 


.004 


2 


,841 


9 
Z 


A Q1 
,4ol 


1 ( 


1 /1 + 
14 


i 


.yys 


9 K91 
Z.OZ1 


q nns 

O.UUo 


3 


.076 


2. 


,908 


9 
Z 


/181 
,4ol 


1 8 + 

lo 


10 


o 

z 


nno 
.uuy 


9 R71 
Z.O I 1 


Q n79 

O.U 1 z 


3 


.143 


2. 


.971 


9 

z. 


coo 
.OZo 


iy 


10 


o 

z 


mo 
.uiy 


9 71 O 

z. / iy 


Q 1 Q9 

o.loZ 


3.207 


3.031 


9 
Z 


c;7Q 


on+ 
zu 


1 7— 
1 l 


z 


D98 
.UZo 


9 7fiQ 
Z. 1 OO 


Q 1 88 

o.loo 


3.266 


3.087 


9 
Z 


.010 


z 


1 


1 


.oUU 


1 7H/1 

1. r y4 


1 707 


1 


.797 


1 


,793 


1 


7n k 

. ( yo 


Q — 

o 


9 + 
Z 


1 

1 


QQQ 
.OOO 


1 Q £ Q 

l.oOo 


1 QSfi 

l.ooO 


1 


.390 


1 


,370 


i 
1 


Q/l /I 
,o44 


A + 
4 


Q — 

o 


1 

1 


97Q 

.Z to 


1 Q91 

l.ozl 


1 Q8Q 
l.OOO 


1 


.391 


1 


,356 


1 


Qni 
oUl 


cr — 



A + 
4 


1 

1 


9 c;o 

.zoy 


1 QQQ 
l.OOO 


1 /1 90 

i.4zy 


1 


.441 


1 


,393 


i 
1 


Qna 
,oUo 








1 

1 


9 c;7 
.ZO < 


i Q7n 
l.o f U 


1 /I 8fi 
1.4oO 


1 


.503 


1 


.442 


i 
1 


Q9Q 

,ozy 


7— 
( 





1 

1 


o c;a 
.zoo 


i /i nf; 
1.4U0 


1 c; I c 
1.040 


1 


.565 


1 


,495 


i 
1 


,o0o 


8 + 
o 


7— 

( 


1 

1 


OKI 
.ZOl 


1 A A Q 
1.44o 


1 fi09 

l.oUz 


1 


.625 


1 


,547 


i 
1 


Q8Q 
.OOO 


o— 

y 


a+ 
o 


1 

1 


.Z04 


1 /i an 
1.4oU 


1 p. KP. 

1.000 


1 


.683 


1 


,598 


i 
1 


A 1 Q 

,41o 


1 n+ 
1U T 


y 


-1 

1 


.ZO f 


1 k 1 £ 

l.olo 


i 7n7 
1. ( U / 


1 


.736 


1 


.646 


1 


A AO 

,44z 


11" 


1 rv-t- 

10 + 


1 


.270 


1.550 


1.754 


1 


.786 


1 


.691 


1 


.471 


12+ 


11- 


1 


.272 


1.582 


1.798 


1 


.833 


1 


.734 


1 


.498 


13" 


12+ 


1 


.275 


1.612 


1.840 


1 


.876 


1 


,774 


1 


.526 


14+ 


13- 


1 


.277 


1.641 


1.878 


1 


.917 


1 


,812 


1 


.552 


15" 


14+ 


1 


.279 


1.669 


1.914 


1 


.955 


1 


,848 


1 


.577 


16+ 


15- 


1 


.281 


1.695 


1.948 


1 


.991 


1 


,881 


1 


.601 


17" 


16+ 


1 


.283 


1.719 


1.979 


2 


.025 


1 


,913 


1 


.624 


18+ 


17" 


1 


.285 


1.742 


2.009 


2 


.056 


1 


,943 


1 


.646 


19" 


18+ 


1 


.287 


1.764 


2.037 


2 


.086 


1 


.971 


1 


.667 


20+ 


19- 


1 


.288 


1.785 


2.064 


2 


.114 


1 


.998 


1 


.688 



20 



Table 5: Experimental B(E1; L -> L - 1) / B(E2; L —> L — 2) ratios (multiplied by 10 5 ) 
[20j of B(E1) and B(E2) values originating from the same level of 226 Th, compared to 
theoretical predictions for O — 45°, 60°, as well as to theoretical predictions by Bizzeti and 
Bizzeti-Sona [20] (labeled as BBS). As in Ref. [201 , asterisks indicate the values used for 
normalization. 





exp 


45° 


60° 


BBS 


8+ 


2.0 (8) 


1.454 


1.457 


1.3 


9" 


1.7 (2) 


1.649 


1.652 


1.3 


10+ 


1.5 (1)* 


1.500* 


1.500* 


1.6* 


11- 


1.7 (1)* 


1.700* 


1.700* 


1.6* 


12+ 


1.6 (1) 


1.544 


1.542 


1.8 


13" 




1.747 


1.746 


1.8 


14+ 


1.4 (1) 


1.585 


1.582 


1.9 


15- 


1.7 (3) 


1.791 


1.789 


2.0 


16+ 




1.622 


1.619 


2.1 


17- 


1.5 (3) 


1.831 


1.829 


2.1 


18+ 




1.656 


1.653 


2.2 


19- 


1.7 (4) 


1.867 


1.865 


2.3 



Table 6: Experimental B(E1; L -> L + 1) / B(E1; L -»• L - 1) ratios |Uj of B(E1) values 
originating from the same level of 226 Ra, compared to theoretical predictions for 0o = 45°, 
56°, as well as to three different theoretical predictions from Ref. jSHl, labelled as R-h, R-I, 
R-II. See section 4 for further discussion. 





exp 


45° 


56° 


R-h 


R-I 


R-II 


1- 


1.85±1.20 


2.116 


2.092 


1.84 


1.98 


1.85 


3~ 


0.87±0.35 


1.451 


1.433 


1.12 


1.31 


0.95 


5- 




1.297 


1.288 








7- 


1.79±1.59 


1.220 


1.215 


0.86 


1.12 


0.99 


9- 


1.27±0.68 


1.172 


1.170 


0.83 


1.10 


1.13 


11- 


1.12±0.79 


1.140 


1.139 


0.83 


1.10 


1.26 


13- 


1.06±0.68 


1.117 


1.117 


0.85 


1.11 


1.35 
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Table 7: Parameter values /3o, m ax where the first derivative of the energy ratios R(L) = 
E(L)/E(2) for the ground state band and the associated negative parity band (s = 1) of 
the Davidson potentials of Eq. (fT6*j) has a maximum, while the second derivative vanishes, 
together with the R(L) ratios obtained at these values (labeled by "var") and the corre- 
sponding ratios of the present model (labeled by "oct"), for several values of the angular 
momentum L. In all cases 0o — 45° has been used. See section 5 for further discussion. 
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Rl 
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1.200 
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0.346 


2+ 
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1.000 
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1.329 


3.030 


3.039 


5" 


1.374 


4.333 


4.351 


6+ 


1.419 


5.797 


5.829 


7" 


1.461 


7.407 


7.459 


8+ 


1.502 


9.154 


9.233 


9~ 


1.541 


11.032 


11.144 
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1.579 


13.034 


13.187 



22 



80- 



60- 



o 



jg 20 
(D 



(a) 

i i 

j \ 

i 

i 


10 R(4) 

R(12) 

R(20) 
R(0 2 ) 


: j 
! I .■ 

j A 
/.•' / 

••/ 1 
-•'/ /' 


i < 
i \ 




/ j 
i 


_ _\ \, 




7 7 









15 30 45 60 75 90 

angle ty (degrees) 



3.0 



2.5 



o 



2.0 



CM 
DO 



1.5 



(b) 



- R{42} 
R{1210} 
R{2018} 

- R{53} 
R{1 31 1} 

■ R{2119} 



15 30 45 60 75 

angle <h (degrees) 



90 



4.0- 



3.5- 



U) 3.0- 
O 



2.5- 



CQ 



2.0- 



1.5- 



R{21} 

R{121 1} 

R{2019} 
R{41} 

R{129} 
R{2017} 



15 30 45 60 75 

angle ty (degrees) 



90 



Figure 1: (color online) (a) Energy ratios R(L) = J^py as a function of the angle O - 
The minima appear as follows: #(4) at 41.19°, #(12) at 42.03°, #(20) at 42.61°, R(0 2 ) 
at 34.99°, R(0 3 ) at 34.13°. (b) Same for the B(E2) ratios #{42} = gg|g^j (40.5°), 

#{1210} = ^ggllff (41.0°), #{2018} = ^gffe^ (42.0°), #{53} = §g|f^j 

(40.5°), #{1311} = B ig2-3-^i-V ( 41 - 5 °)> ^{2119} = ^e^Z 1 !-) (42.0°). After each 
ratio, the position of the maximum (with accuracy ±0.5°) appears in parentheses, (c) 
Same for the B(E1) ratios #{21} = fjjgg^j (40.0°), #{1211} = (40.5°), 

#{2019} = ^fftlff (41-5°), and for the #(#3) ratios #{41} = g^lfe^ (41-0° 

^{129} = §ft|^ (41-0°), #{2017} = ^ggt^y 



B(£3;3-^0+) 

(41.5°). Again, each ratio is 
followed by the position of the maximun (with accuracy ±0.5°) in parentheses. 
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angular momentum L angular momentum L 

Figure 2: (color online) (a) Experimental energy ratios R(L) = E(L)/E(2f) for 220 Th 
HH, 222 Th |2H], 224 Th j2H], 226 Th |3Dj, 228 Th jHU, 230 Th |3j, 232 Th (S21IS2!, and 234 Th 
[32j, compared to theoretical predictions for <fi = 45° and <ft — 60°. (b) Same for 218 Ra 
jSlEni, 220 Ra j2Zl, 222 Ra [321 EE!, 224 Ra |S2lEni, 226 Ra |32l Eg, and 22§Ra E2j, compared 
to theoretical predictions for = 45° and = 56°. 
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Figure 3: (color online) (a) Experimental energy levels of the ground state bands of 220 Th 

|27j, 222 Th [2E], 224 Th ESI, 226 Th Em, 228 Th eu, 230 Th pa, 232 Th E3E31, and 234 Th E2l, 

as a function of the neutron number, (b) Same as (a), but for the associated negative parity 
bands, (c) Same as (a), but for 218 Ra [MJ ESI, 22 ° Ra EZj, 222 Ra [S3 ES, 224 Ra El EE], 
226 Ra E21EH1, 228 Ra [32], and 230 Ra [32] ■ (d) Same as (c), but for the associated negative 
parity bands. 
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Figure 4: (color online) (a) Experimental energy staggering AE (Eq. flHBJ)), normalized to 
the 2f state of each nucleus, for 220 Th [27], 222 Th j2H], 224 Th [2H], 226 Th [30 , 228 Th f^H , 
230 Th |32l, 232 Th (221ES1, and 234 Th [32], as a function of the angular momentum, (b) Same 
as (a), but for 218 Ra HUES], 220 Ra |JJ7J, 222 Ra E21IM], 224 Ra [32lll>], 226 Ra [32 EI], 228 Ra 
]S2], and 230 Ra ]S2]. 
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Figure 5: (color online) (a) Experimental B(E1;L — > L — 1) / B(E2;L — > L — 2) ratios 
(multiplied by 10 5 ) [20] of B(E1) and B(E2) values originating from the same level of 226 Th, 
compared to theoretical predictions by Bizzeti and Bizzeti-Sona j20] (labeled as BBS), as 
well as to predictions of the present work for </>n = 45°, 60°. The ratios corresponding 
to L = 10 and 11 have been used for normalization, as in Ref. [20J. (b) Experimental 
B(E1; L -> L + 1) / B(E1; L -> L - 1) ratios |U] of values originating from the 

same level of 226 Ra, compared to three different theoretical predictions from Ref. jSH] 
(labeled as R-h, R-I, R-II), as well as to theoretical predictions of the present work for 
0o = 45°, 56°. See section 4 for further discussion. 
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Figure 6: (color online) Energy ratios R(L) = E(L)/E(2) for the ground state band and 
the associated negative parity band (s = 1) of the Davidson potentials of Eq. (fTBJ). selected 
through the variational procedure of section 5 (labeled by "var"), compared to the ratios 
provided by the present work (labeled as "oct"), as a function of angular momentum. The 
limiting cases corresponding to Davidson potentials with @ Q = (labeled as "/3n = 0" ) and 
/3q —> oo (labeled as "rot") are also shown for comparison. In all cases, 0o = 45° has been 
used. See section 5 for further discussion. 
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